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The point defect thermodynamics in a general family of binary compounds, including B2 com¬ 
pounds as a specific representative, are classified by way of two non-trivial energy parameters. The 
scheme is applied to published ab initio defect formation energies, and the variety of resulting phe¬ 
nomena is demonstrated. Further, by introducing model assumptions the consequences for the active 
diffusion mechanisms are deduced. It is shown that particularly for the off-stoichiometric case, the 
assumed prevalence of either the six-jump cycle or the triple-defect mechanism has to be reconsid¬ 
ered, as a number of qualitatively different mechanisms emerge as likely candidates for the dominant 
effect. Two of those, the 4-|-2-jump cycles and the waltzing-step mechanism, are introduced here. 


I. INTRODUCTION 

In pure elemental crystalline systems, the question of 
the thermal point defects’ nature as well as of the ac¬ 
tive mechanisms for self diffusion can be considered as 
largely solved. Specifically for metals, the native point 
defect is the vacancy, whose random, thermally activated 
walk is also responsible for the stochastic translation of 
the atoms leading to self diffusion, while for semiconduc¬ 
tors also diffusion via self-interstitials is plausible. For 
ordered compounds, the situation is much more com¬ 
plicated. Neglecting interstitials, already for a binary 
compound composed of two sublattices four fundamental 
point defects are possible, namely vacancies and antisite 
atoms on either sublattice. Further, the concentrations of 
these defects for a given system as a function of temper¬ 
ature are not independent variables, as experimentally 
the composition of the system has to be considered as 
fixed. Several theories to solve this issue have been put 
forward.lilf^ 

The mechanisms responsible for diffusion in an ordered 
compound are still controversial. An uncorrelated ran¬ 
dom walk of a vacancy composed of nearest-neighbour 
jumps, as familiar from elemental metals, can be ruled 
out, as it would destroy the ordered structure during its 
migration. Specifically the B2-structure and analogous 
compounds, where a given atom’s nearest neighbours are 
exclusively on the sublattice occupied by the other el¬ 
ement, a vacancy could only move by jumps to farther 
neighbours within one sublattice without introducing an¬ 
tisite atoms. While such next-nearest neighbour jumps 
are a priori entirely possible, it was recognized early on 
that they would correspond to large migration energies. 
H. B. Huntington is attributecP the first suggestion of a 
mechanism where a vacancy performs a correlated walk, 
so that it annihilates at later steps the antisite atoms cre¬ 
ated during earlier steps, leading to diffusion via nearest- 
neighbour jumps while conserving the atomic order in 
the long run. Apart from next-nearest neighbour jumps, 
these so-called six-jump cycles and the later proposal 
of the triple-defect mechanisnP are the candidates that 
are conventionally considered for explaining experimen¬ 


tal diffusion datcP or for predicting the active mechanism 
from theoretical calculation^. However, it is conceivable 
that other mechanisms have hitherto been overlooked, 
which is sometimes conceded explicitly (e.g.. Ref. El). 

The aim of this article is to explore the range of qual¬ 
itatively different behaviours with respect to constitu¬ 
tional and thermal defect structures and, related to this 
issue, the distinct conceivable diffusion mechanisms, with 
special consideration to the B2 structure as one of the 
simplest intermetallic structures. It will recapitulate the 
results of a grand-canonical theory of defect concentra¬ 
tions, propose a re-parametrization of the effective de¬ 
fect formation energies as reported by theoretical calcu¬ 
lations, derive the phase diagram of dominant constitu¬ 
tional and thermal defects with respect to these energy 
parameters, and classify published ab-initio defect ener¬ 
gies for a variety of B2 systems accordingly. Further, an 
enumeration of the plausible diffusion mechanisms me¬ 
diated by nearest-neighbour vacancy migration will be 
given, along with an evaluation of their respective en¬ 
ergies for above parametrization under two simple mod¬ 
els. This analysis will yield the surprising result that, in 
contrast to prevailing assurance, neither six-jump cycles 
nor the triple-defect mechanism are expected to be re¬ 
sponsible for diffusion in many prominent systems, but 
one of two qualitatively different mechanisms proposed 
here. Finally, extant knowledge on diffusion in NiAl and 
CoGa will be discussed in the light of the insights out¬ 
lined above. 


II. POINT DEFECT THERMODYNAMICS 
A. Grand-canonical theory 

Here the main results of Fahnle and coworkers’ grand- 
canonical theory for the point defect concentrations in 
ordered structures will be recapitulated (for a thorough 
account see, e.g.. Refs. [3] and [5] and the references cited 
therein). It deals with the point defect concentrations 
in ordered structures neglecting defect interactions, in 
its initial formulation for an ordered binary structure of 
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AB stoichiometry with two sublattices. Strictly speak¬ 
ing, the latter condition can be relaxed to require only 
chemical, but not crystallographic equivalence between 
one constituent’s native sites. This encompasses most of 
Strukturbericht type B structures such as B1 (NaCl), B2 
(CsCl), B3 (ZnS), or LIq (CuAu), to name only the most 
prominent, along with B20 (FeSi) as an example with 
eight sublattices, which, however, correspond to only two 
chemically inequivalent sites. While quantitative gener¬ 
alizations for situations with higher complexity, such for 
as DO 3 , have been giverfi^, the qualitative effects can di¬ 
rectly be derived from the special case treated here. Note 
also that while there exist generalizations to treat in¬ 
teractions between defects in an approximate wajff^, ne¬ 
glecting interactions is exact in the limit of small defect 
concentrations as long as no complex defects appear and 
the system is thermodynamically stable (i.e., no phase 
separation by defect agglomeration occurs). 

Four basic point defects are considered, to wit vacan¬ 
cies or anti-structure atoms on either sublattice (intersti¬ 
tials are neglected for simplicity, as their equilibrium con¬ 
centration has been recognized to be negligible in com¬ 
parison to vacancies and antisites in intermetallic com¬ 
pounds). The salient point of the theory is that in or¬ 
dered compounds the notion of a given basic point de¬ 
fect’s formation energy is ill-defined, as due to particle 
number conservation any such defect has to be generated 
in conjunction with at least another one. Therefore, one 
chooses a grand-canonical setting and considers so-called 
grand-canonical defects, where the particle numbers are 
not conserved and whose formation energies are defined 
up to one degree of freedom that has the role of a chem¬ 
ical potential. For any given temperature this chemical 
potential is implicitly defined so as to give the correct 
composition. 

The main result is that in the non-degenerate low- 
temperature limit a single type of defect dominates in 
the off-stoichiometric case (termed the constitutional de¬ 
fect). Therefore the chemical potential converges to some 
specific value, and the equations describing the concen¬ 
tration of the various basic defects simplify to Arrhenius- 
type expressions defined by so-called effective formation 
energies, which do not depend on temperature or com¬ 
position. However, the pre-factor of these expressions 
includes a concentration-dependent term in addition to 
the conventional electronic and vibrational contributions 
to defect formation entropies. Analogously, in the case of 
perfect stoichiometry two basic point defects appear in 
coupled concentrations as the dominant low-temperature 
thermal defects, and the same simplifications ensue. For 
compositions and temperatures where the numbers of 
constitutional and thermal point defects are compara¬ 
ble, the full set of coupled non-linear equations has to be 
solved. In the following it is assumed that either consti¬ 
tutional or thermal defect concentrations dominate over 
the other. 

To define the nomenclature, let A and B be the con¬ 
stituent elements, with their native sublattices a and fi. 


Naturally, an a vacancy is an empty site on the a sublat¬ 
tice, and an a antisite is a site on the a sublattice occu¬ 
pied by a B atom, which will be called B anti-structure 
atom. To motivate the behaviour reported above, de¬ 
fine e^+fj. a.s the grand-canonical energy difference when 
removing one A atom from its native a-sublattice, anal¬ 
ogously — p, for removing B from /3, further e® -|- 2/i 
and e-p — 2 /i as the energy differences for replacing one 
A atom on its native sublattice by B and vice versa, re¬ 
spectively. Note that when calculating a given configu¬ 
ration’s total energy, the term proportional to /i indeed 
depends only on the composition, so that introducing 
any number of defects while keeping the particle number 
constant is independent of /i, the chemical potential for 
exchanging B for A. Now choose /i so that each of above 
energy differences is positive (this is possible, as other¬ 
wise the ordered structure would be thermodynamically 
unstable). If now temperature is lowered for fixed /x, the 
composition of the system evolves towards stoichiometry 
(each defect gets progressively less probable). To keep 
the system at a certain off-stoichiometric composition, /r 
has to be adjusted in step with temperature (clearly the 
composition varies monotonically with /i, so there is a 
unique solution). It is now obvious that for a certain off- 
stoichiometric composition there has to be a defect whose 
energy goes to zero as temperature goes to zero. This de¬ 
fect is the constitutional defect with effective formation 
energy equal to zero, and the other defects’ effective for¬ 
mation energies follow from the corresponding //. In the 
case of stoichiometry on the other hand, for a general 
choice of /i one defect will have the smallest formation 
energy, so its concentration will dominate over the other 
defects for small temperature, which violates the condi¬ 
tion of stoichiometry. As a consequence, in this case /i 
has to be chosen so that the minimal effective formation 
energy is shared by two complementary defects (i.e., ei¬ 
ther both vacancies or both antisites, or a vacancies and 
/3 antisites or vice versa). 

B. Defect formation energy parametrization and 
phase diagrams of defect types 

This section proposes a reparametrization of the de¬ 
fect formation energies in order to explore the range of 
resulting defect structures. To this end define 

^0 = £« + + (e;a + ea)/2 (la) 

C = ((^a + e^) — (e^ + ea)/2)/£^o (lb) 

’I = ((Ca ~ e^) + ~ ea)/ 2 )/l?o- (Ic) 

Note that for computing these quantities, the grand- 

canonical energy differences for any choice of chemical 
potential /i can be used, as the corresponding contribu¬ 
tion cancels in each case. Therefore, these three param¬ 
eters define the defect energetics of a given system in a 
unique and non-redundant way. The positive parameter 
Eq has the dimension of an energy and its only action 
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FIG. 1. Phase diagrams of constitutional and dominant thermal point defects of an AB-system for stoichiometry (left) and 
A-excess (right), where the dashed line separates the regions of constituional antisites (above diagonal) and vacancies (below). 
Dominating thermal excitations are coded by colour, where orange and gray are vacancy and antisite pairs, blue and yellow 
are A and B triple defects, and green and pink are antisite and vacancy annihilation, respectively. 


corresponds to a rescaling of temperature, while ^ and 
1 ] are dimensionless and vary from —1 to 1. In the fol¬ 
lowing, the constitutional and thermal defect structures 
correspond to a given choice of (^, rj) will be explored. 

For the stoichiometric case, the four possible thermal 
excitations of the site occupations are (1) a vacancy pair, 
wherein the crystal is enlarged by one cell, generating one 
a and one /3 vacancy, (2) an antisite pair, wherein an A 
atom is swapped with a B atom, generating one a and 
one P antisite, (3) a triple defect in A, wherein the crystal 
is enlarged by one cell and an A atom is moved into 
the resulting /3 vacancy, generating two a vacancies and 
one /3 antisite, and lastly (4) analogously a triple defect 
in B. In the case of vanishing defect concentrations, the 
configurational entropy corresponding to one of above 
excitations is equal to the number of affected sites times 
kn log(iV), where N is the number of sites per sublattice 
(this constant factor will be suppressed in the following). 
The formation energies and entropies are tabulated in 
Tab. The dominant thermal excitation resulting for a 
specific (^, T]) is the one with the smallest cost in energy 
per gain in entropy, as given in the last column of the 
table. The corresponding phase diagram is drawn in Fig. 
[^(left), with the single critical point being at (1/3,0). 

For off-stoichiometry only the case of A-excess shall be 
treated, as the results for B-excess follow by symmetry. 
Here the primary distinction is the nature of the consti¬ 
tutional defects. An excess of two A atoms can be ac¬ 
commodated either by two /3 vacancies or one j3 antisite. 
Therefore, if 2e^ — = i?o(C — ??) > 0 the constitutional 


excitation 

energy E entropy S E/SEq 

vacancy pair 

€ck “T 

2 (1 + 0/4 

antisite pair 

-1- 
E/3 + 

2 (1-0/2 

A triple defect 

2ea -t eja 

3 (1 + 0/3 

B triple defect 

+ e® 

3 (1-0/3 


TABLE I. Formation energies, entropies and formation energy 
per entropy increase of thermal occupational excitations for 
stoichiometry 


defect will be the (3 antisite, otherwise the j3 vacancy. 
The question of the thermal excitations is analogous to 
the stoichiometric case with the following complications: 
First, in addition to above four excitations there exists 
in the case of constitutional antisites the possibility of 
(5) antisite annihilation, wherein the crystal is enlarged 
by one cell and one constitutional A anti-structure atom 
is moved into the a vacancy, generating two j3 vacan¬ 
cies and annihilating one /3 antisite, and in the case of 
constitutional vacancies analogously the possibility of (6) 
vacancy annihilation, wherein one A atom moves into a 
constitutional /3 vacancy and the crystal is shrunk by 
one cell. These additional excitations have been called 
interbranch defects by Korzhavyi et al.E^. Second, as 
the thermal defect concentration is by assumption much 
smaller than the constitutional defect concentration, the 
thermal excitation (or annihilation) of point defects that 
exist already as constitutional defects does not affect the 
entropy of the system. The relevant expressions are given 












4 


excitation 

energy E entropy S 

E/SEo 

vacancy pair 

€a -T 

2 

(1 + 0/4 

antisite pair 


1 

1-C 

A triple defect 

SEq + 

2 

(1 + 0/2 

B triple defect 

+ Ca 

3 

(1-0/3 

antisite annihil. 


2 

(C - 0/2 


excitation 

energy E entropy S E/SEq 

vacancy pair 

H" 

1 (1 + 0/2 

antisite pair 

4 + e® 

2 (1-0/2 

A triple defect 

2ea + 

3 (1 + 0/3 

B triple defect 

2e5 + 

1 1 — 77 

vacancy annihil. 

— 2e5 + 

1 -C + 77 


TABLE II. Formation energies, entropies and formation en¬ 
ergy per entropy increase of thermal occupational excitations 
for A-excess, for the cases of constitutional a antisites (above) 
or p vacancies (below). 


in Tab. [IT] Note that from the argumentation in Sect. |II A| 
it follows that the dominant thermal point defect is the 
one with the lowest effective formation energy after the 
constitutional point defect with vanishing effective for¬ 
mation energy. Therefore only those thermal excitations 
that introduce just one kind of thermal point defect in 
addition to the constitutional defect can be expected to 
dominate. This is in agreement with the dominant ex¬ 
citations that follow from Tab. |T^ For instance, in the 
case of constitutional antisites E/S for the vacancy pair 
is equal to the mean values of E/S for A triple defects 
and antisite annihilation, so one of those has smaller E/S 
and vacancy pairs will hence not occur for any choice of 
(^, 77 ). The resulting phase diagram is drawn in Fig. [b 
( right). Here the critical points are given by (3/5, — 1/^ 
and ( 0 , 1 / 2 ). 

To demonstrate that the theoretical richness of defect 
thermodynamics is actually realized, published values of 
defect energetics in a range of B2 systems obtained by 
ab-initio density functional theory calculations have been 
compiled and are reported in Tab. m The qualities of 
the respective calculations differ, which is evidenced for 
instance by the scatter of the values corresponding to 
NiAl, a popular system for defect thermodynamics in¬ 
vestigations. However, plotting the energy parameters in 
the defect structure phase diagrams of Fig. shows that 
the differences between different systems are larger than 
these deviations, so that the qualitative conclusions seem 
sound. Note that the resulting classification is as exact 
as the input from the ab-initio calculations. 

The two panels of Fig. show that indeed all prin¬ 
cipally possible defect structures are realized for some 
(^, 77 ): For the stoichiometric case the thermal excita¬ 
tion has to be composed of one point defect that in itself 
would lead to A-excess and one to B-excess. Of each 
there are two possibilities, so there result four distinct 
thermal excitations. On the other hand, the case of off- 
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Eo 

e 

V 

Ref. 


(eV) 

(eV) 

(eV) 

(eV) 

(eV) 




FeAl 

0.81 

4.01 

0.67 

0.67 

5.49 

0.756 

-0.583 

M 


1.06 

3.46 

0.99 

0.99 

5.51 

0.641 

-0.436 

M 

RuAl 

1.46 

3.07 

1.03 

1.03 

5.56 

0.629 

-0.290 

ESI 


1.45 

3.21 

1.45 

2.02 

6.39 

0.458 

-0.321 

m 

CoAl 

1.29 

2.23 

1.29 

2.91 

5.62 

0.253 

-0.311 

M 

NiAl 

0.94 

2.13 

0.94 

2.21 

4.64 

0.322 

-0.393 

ESI 


0.74 

1.97 

0.74 

2.36 

4.26 

0.272 

-0.479 

EH 


0.79 

1.74 

0.79 

2.85 

4.35 

0.163 

-0.455 

EH 


0.55 

1.53 

0.55 

2.08 

3.40 

0.225 

-0.514 

EZI 


0.73 

1.62 

0.73 

2.68 

4.05 

0.159 

-0.461 

ESI 

PdAl 

1.22 

1.58 

1.22 

2.31 

4.56 

0.227 

-0.198 

m 

Pdin 

0.31 

0.97 

0.31 

1.81 

2.34 

0.094 

-0.603 

EH 

NiTi 

1.09 

1.74 

0.65 

0.65 

3.48 

0.626 

-0.187 

ESI 

FeCo 

1.65 

1.68 

0.15 

0.15 

3.48 

0.914 

-0.009 

ED] 


1.78 

1.78 

0.17 

0.17 

3.73 

0.909 

0.000 

M 


TABLE III. Effective defect formation energies for stoichiome¬ 
try in AB nomenclature as taken or computed from published 
ab-initio values, along with corresponding energy parameters. 


stoichiometry corresponds to one point defect having zero 
effective formation energy, so treating without loss of gen¬ 
erality only A-excess, two choices are possible. For each 
of these, one of the remaining three point defects will 
have the lowest effective formation energy and therefore 
dominate, giving six possible cases of defect structures. 
What is remarkable is that, with the exception of va¬ 
cancy pairs, all these options are realized already by the 
small subset of intermetallic systems in Tab. |HI| For com¬ 
pounds of increasingly ionic character, vacancies will de¬ 
crease in formation energy compared to antisite defects, 
which would favour thermal vacancy pairs. In fact, for 
ionic compounds with comparative sizes of cations and 
anions, such as for the prototypical CsCl, the Schottky 
defect (i.e., vacancy pairs) is known to be the dominant 
thermal excitatioiP^. 

It is customary to classify the thermodynamics of AB 
compounds into antistructure-type systems and triple¬ 
defect systems. This is due to NeumanrP^, who observed 
the dominant defect type to be connected to the forma¬ 
tion enthalpy of the compound. Hybrid behaviour was 
predicted for borderline cases such as FeAl. In further 
work, the distinction was interpreted as the opposite ex¬ 
tremes of parameter choices in general thermodynamical 
models (e.g. Ref. [M|) . with the label “hybrid behaviour” 
assigned to intermediate cases. In most cases, it is un¬ 
derstood that the distinction refers to the nature of con¬ 
stitutional defects, where “anti-structure system” means 
that deviations in both directions from stoichiometry are 
accomodated by constitutional antisites, while in “triple¬ 
defect systems” only, say, A-excess is accomodated by A 
anti-structure atoms, while A-deficiency is accomodated 
by a vacancies. However, sometimes it is also implied 
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that, e.g., in anti-structure systems not only the consti¬ 
tutional defects, but also the “governing” ther mal d efects 
are antisites (equally for triple-defect systemswhile 
in other cases the distinction is assumed to be defined 
directly in terms of the dominant thermal excitations in¬ 
stead of the way off-stoichiometry is accomodatecpSl. 

Fig.[ 2 shows that care has to be taken, as above defini¬ 
tions are by no means equivalent or sometimes even well- 
defined: First, the dominant thermal excitation depends 
on whether the composition is A-deficient, stoichiomet¬ 
ric, or A-excessive. Therefore, there are two different 
consistent possibilities to interpret a given system being 
classified as a triple-defect system in terms of the domi¬ 
nant thermal excitation, namely either only for the stoi¬ 
chiometric case, corresponding to \ri\ > max( 6 ^— 2 , - 3 ^- 1 - 
l)/4, or for all cases, corresponding to \r]\ > (1 -f |^|)/2. 
In contrast, the common understanding of constitutional 
antisites on one side of stoichiometry and constitutional 
vacancies on the other would correspond to \r]\ > |^|. 
Such a distinction in terms of constitutional defects seems 
most natural, so that calling this situation the triple¬ 
defect case is unfortunate, as a triple defect is always 
a thermal excitation. It is conceivable that this is¬ 
sue is partly responsible for the inconsistent interpreta¬ 
tions. Also, claims of hybrid behaviour, e.g., for FeAl 
can have resulted from triple defects being the dominant 
thermal excitation for Fe-rich compositions, while off- 
stoichiometry is accomodated by antisites, so that the 
most abundant thermal point defect is the Fe vacancy. 
Concludingly, it seems most expedient to use a top-level 
classification of AB compounds according to the con¬ 
stitutional defects of anti-structure accomodation with 
^ > \r]\, mixed accomodation (anti-structure atoms on 
one side of the stoichiometry, vacancies on the other) for 
\ri\ > 1^1 and vacancy accomodation for the hypothetical 
case of ^ < — 1 ? 7 |, which apparently has not been reported 
yet. 

The effective defect formation energies are given by 


£«/ — 1/4 -f C/4 + 2r]/5 -I- /i (2a) 

e}/Eo = 1/4 + C/4 - 277/5 - 77 (2b) 

el/Eo = 1/2 - C/2 - 77/5 + (2c) 

e^/Eo = 1/2 - C/2 + 77/5 - 2fi (2d) 

with 

vacancy pairs: fi = — 277/5 (3a) 

antisite pairs: fi = rj/10 (3b) 

A triple defects: /r = (1/4 — 3C/4 — ?7/5)/3 (3c) 


B triple defects: /z = (—1/4-|-3C/4 — 77 / 5)/3 (3d) 

for the regions of different dominant thermal excitations 
in stoichiometry and 

structural antisites: 77 = 1/4 — C/ 4 -l-? 7/10 (3e) 
structural vacancies: 77 = 1/4 -|- C/4 — 277/5 (3f) 

for A-rich off-stoichiometry. Note that always the mirror 
case for a given expression (i.e., exchanging the roles of 


A and B as well as those of a and /?) can be obtained 
by flipping the signs of 77 and 77 . It is easy to verify that 
Eqs. are consistent with Eqs. ii). and that Eqs. @ in¬ 
deed lead to vanishing effective formation energies of the 
respective constitutional defects in the off-stoichiometric 
case and to equal effective formation energies of the prim¬ 
itive defects making up the dominant thermal excitations 
in the stoichiometric case. 

It is worth noting that for a given system not neces¬ 
sarily all primitive defects actually exist. For instance, 
in NiAl the A1 vacancy is separated only by a barrier 
of 0.35 eV against the decay in to a N i vacancy and a Ni 
anti-structure atom (see Sect. IVA). For higher asym¬ 
metries, the barrier could vanish, so that the isolated 
vacancy would exist only as an extremely anharmonic os¬ 
cillation on phonon timescales. However, for the present 
purposes it would be correct to use the energy of the 
nearest-neighbour vacancy-antisite pair in place of the 
unstable vacancy, with an appropriate interpretation of 
the resulting point defect concentrations. 

Finally, for most systems of interest the criteria for 
the stoichiometric case, i.e., that thermal defects are 
much more abundant than constitutional defects, are 
quite stringent. Take for instance NiAl with an effec¬ 
tive formation energy of thermal antisites around 0.75 eV 
according to Tab. m For an intermetallic compound de¬ 
viations of at least 0 . 1 % from nominal stoichiometry have 
to be realistically expected in the experiment. Only at 
the congruent melting temperature of 1911K the ther¬ 
mal defect concentration reaches 1 %, so that the stoichio¬ 
metric limit can be considered as valid. In contrast, the 
off-stoichiometric limit will be appropriate for the major¬ 
ity of experimental cases in well-ordered systems. Only 
for weakly-ordered systems such as FeCo, which display 
an order-disorder transition within the solid phase, large 
concentrations of thermal defects will exist at intermedi¬ 
ate temperatures. 


III. DIFFUSION 
A. Transition state theory 

Experiment has proven that the hopping rate, and as a 
consequence the diffusion constant, of a single jump dif¬ 
fusion mechanism such as vacancy diffusion in elemental 
crystals can often be described well by an Arrhenius-type 
expression. Indeed, transition state theory as formulated 
by VineyarcP^ leads to exactly such an expression for the 
hopping rate, where the activation energy is equal to the 
difference between saddle point energy and initial energy 
in the 37V-dimensional energy landscape of the A^-atom 
subsystem considered explicitly, and the frequency pref¬ 
actor is the product of all 3N eigenfrequencies in the 
initial stable state divided by the 3N — 1 real eigenfre¬ 
quencies in the saddle point, provided the harmonic ap¬ 
proximation is valid in both states. In addition to this 
condition, whose failure can result in deviations from Ar- 
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rhenius behaviour, transition state theory assumes that 
once a system crosses over the saddle point, it continues 
until the neighbouring local minimum, where it thermal- 
izes. If the time scales of the jump duration and the 
mean residence time in a local minimum differ enough 
and the saddle points in SiV-dimensional space are nar¬ 
row enough so that an immediate further jump is im¬ 
probable, then N can be chosen so that the subsystem 
behaves microcanonically during the jump, i.e., does not 
reverse, while the coupling to the heat bath constituted 
by the other atoms leads to thermalization in the final 
local minimum. These conditions are normally fulhlled 
for bulk vacancy diffusion. 

As discussed below, in compounds often correlated 
jump sequences are relevant for diffusion. There, the 
system has to pass through a number of local energy 
minima, called stable states from now on, for complet¬ 
ing the jump sequence. In each of these it will rest long 
enough for thermalizing, and therefore can potentially 
pass back over the energy apex along the jump path. As 
a consequence, the above account has to be extended for 
treating diffusion in compounds. 

Here the comparatively simple case of diffusion via a 
unique string of stable states will be considered. Let 
the initial stable state be the energy minimum along the 
string, called the ground state, and denote its energy by 
E^. It is separated by a saddle point with energy Ef 
from the next stable state with energy A™, up to the 
final stable state M, with the corresponding products 
of eigenfrequencies w™ and wf. According to transition 
state theory the flux between states i and z -I- 1 is given 
by 










i±L 


(4) 


where q is the concentration in the respective stable 
states. For obtaining the steady-state flux identify all 
fluxes with <&, solve Eq. Q for a+i and perform repeated 
substitution, yielding 

^ CotU^e-^o/kBT _ 

Comparing Eqs. Q and ([^ shows that the simple Ar¬ 
rhenius expression for diffusion via direct jumps between 
stable states has to be replaced by the reciprocal sum 
of such expressions for diffusion via intermediate stable 
states. However, for realistic choices of wf and Ef the 
errors introduced by assuming a simple Arrhenius be¬ 
haviour with the highest saddle point energy as activation 
energy are minute, specifically it is exact in the small- 
temperature limit. This qualitative conclusion agrees 
with detailed calculations of, e.g., the six-jump cycles in 
B2-compounds, whe re the re are several paths connecting 
initial and final stat^SHMl^ 


B. Diffusion mechanisms 


Arguably the simplest class of ordered compounds is 
constituted by the binary compounds of AB stoichiom¬ 
etry where all nearest-neighbour sites of an a site are /3 
sites and vice versa. Here a list of conceivable diffusion 
mechanisms for these compounds will be given. In Sect. 
|HI C| it will be argued that under plausible assumptions 
this list is exhaustive. A two-dimensional square lattice 
with alternating occupations will be used as a model for 
illustrating the mechanism. For actual intermetallic com¬ 
pounds, the B2 structure is by far the most prominent 
representative of this class of compounds, and special 
consideration will be devoted to it below. 

Historically, with the growing acceptance of stochas¬ 
tic nearest-neighbour vacancy jump diffusion being the 
dominant mechanism of diffusion in elemental metals, 
it was recognized that such a mechanism is incompat¬ 
ible with the observed stability of the ordered struc¬ 
ture during diffusiorP^. Apart from the possibility of 
any element performing jumps directly within its na¬ 
tive sublattice, corresponding to jumps to second-nearest 
neighbour sites, correlated jump mechanisms composed 
of nearest-neighbour jumps have been proposed to obvi¬ 
ate the high migration energies assumed for these farther 
jumps. These correlated jump sequences are the six-jump 
cyclePl and the divacancy mechanism as believed to be 
active in ionic compounds, where it suggests itself al¬ 
ready because of its favourable Coulomb energy. In the 
case of intermetallic compounds, the divacancy mecha¬ 
nism is often called triple-defect mechanisnP, as for large 
thermodynamical asymmetries a localized triple defect as 
opposed to the divacancy can be the ground state for this 
jump mechanism. At large deviations from stoichiome¬ 
try diffusion via percolating antisite defects can become 
possible, termed the anti-structure bridge mechanisrrPi^. 


The predominant opinion today is that diffusion will 
happen via the anti-structure bridge mechanism if possi¬ 
ble, otherwise via one of the correlated jump mechanisms, 
i.e., the six-jump cycle or the triple-defect mechanisnpZl. 
While according to calculations the migration energies for 
jumps between second-nearest neighbours are not as high 
as to be ruled out a priori (compare, e.g., Ref. IB5lfor NiAl, 
which is arguably the most prominent system for these 
kinds of studies), the constituents’ respective activation 
energies of diffusion often seem to be coupled, which hints 
towards correlated jump mechanisms. The following sys¬ 
tematic discussion will be restricted to jumps between 
nearest-neighbour sites at low temperatures so that no 
thermal disorder has to be considered and differentiating 
between clear-cut diffusion mechanisms is possible. 

The first qualitative distinction is whether the mecha¬ 
nism can be active in an ideally ordered lattice or whether 
it needs percolating constitutional defects. For the lat¬ 
ter possibility in the case of constitutional antisites this 
is the anti-structure bridge mechanism, where a vacancy 
jumps from the majority sublattice onto a constitutional 
antisite and back onto the majority sublattice. It is il- 
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FIG. 2. Diffusion on percolating antisite cluster — anti¬ 
structure bridge mechanism, a sites have diamond shape, 
/3 sites have square shape; A atoms are blue, B atoms red, 
vacant sites are white. The thermal vacancy’s movement via 
a sites and constitutional /3 antisites leads to diffusion of the 
majority A component. 

lustrated in Fig. for A-excess (in the following called 

anti-structnre bridge mechanism, due to the j3 anti¬ 
sites acting as bridges). Evidently no additional disor¬ 
der is introduced, therefore this mechanism is commonly 
assumed to have a small activation energy and be the 
dominant mechanism responsible for diffusion of the ma¬ 
jority component once percolation permits. Ref. [34] re¬ 
ports it to set in for the B2 structure at a composition of 
around 54.9(2)% of A atoms (i.e., where 10% of /? sites 
are occupied by A atoms). A more precise value for the 
percolation threshold in the academic case of ideal ran¬ 
dom occupancies corresponding to no defect interactions 
is 54.882225(5)9pl. 

Also for the mirror situation of constitutional vacancies 
a low-energy percolation mechanism for majority species 
diffusion exists, which is illustrated in Fig. For the 
case of constitutional /3 vacancies, the diffusing defect is 
a thermal /3 antisite. The elementary step here consists in 
a nearest neighbour of the anti-structure atom jumping 
into a constitutional vacancy, followed by the initial A 
anti-structure atom moving into the transient a vacancy. 
In effect, the thermal antisite defect migrates on the con¬ 
stitutional /3 vacancies, just as in the case of the anti¬ 
structure bridge mechanism a thermal vacancy diffuses 
via constitutional antisites. The connectivity conditions 
are the same, so the same percolation threshold of 10% 
vacancy concentration on the minority sublattice applies 
for the B2 structure. This corresponds to a deviation 
from stoichiometric composition of only 2.5%, half of the 
value necessary for the anti-structure bridge mechanism. 
This mechanism has probably first been considered by Xu 
and van der Ven, who found it to be active in their sim¬ 
ulations of Al-rich NiAP. They called it anti-structure 
bridge sequence for A1 atoms, which is unfortunate as 
it does not allow to distinguish this mechanism from a 
proper anti-structure bridge mechanism on thermal an¬ 
tisites as proposed by Herzig and Divinski for the same 
systeirP^. Therefore, here it shall be termed vacancy 
bridge mechanism, as the constitutional vacancies pro¬ 
vide the bridges for the diffusing thermal antisite defect, 
just as the constitutional antisites provide the bridges 
for the diffusing thermal vacancy in the anti-structure 
bridge mechanism. Note that the vacancy bridge mech¬ 
anism can also be considered a variant of the classical 
interstitialcy mechanism of diffusiorP^, where the perco- 


FIG. 3. Diffusion via percolating vacancy cluster — vacancy 
bridge mechanism. Interpretation of symbols as in Fig.The 
movement of the thermal P antisite between the constitutional 
P vacancies leads to diffusion of the majority A component. 
The defect configuration and migration path mirror those of 
Fig.HI to demonstrate the analogy. 

O u O ^ OoO O u O O m O 
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FIG. 4. Gorrelated propagation of neutral defect — divacancy 
mechanism. 

lating vacant sites take the role of the interstitial sites’ 
lattice. 

For discussing diffusion mechanisms where a well- 
defined defect migrates in an otherwise ideal environ¬ 
ment, it is useful to distinguish by way of the diffus¬ 
ing defect’s charge, as this quantity is conserved over all 
stages of the mechanism. To this end assign the charge 
of -|-1 to an A atom irrespective of whether it is on its 
native sublattice or not, 0 to a vacancy, and — 1 to a B 
atom. In other words, the charge as defined here is the 
difference of the numbers of A and B atoms counted over 
the lattice cells that cover the migrating defect. 

For each defect charge, there exists a qualitatively 
unique string of states, i.e., assignments of atoms to lat¬ 
tice sites, linked by nearest-neighbour vacancy jumps, 
that generates the smallest number of antisite defects and 
that therefore is expected to correspond to the smallest 
migration energy. However, depending on the system, 
its actual realization can take a number of forms, dis¬ 
tinguished by the geometry of the jump path, or which 
states along the path correspond to the ground and apex 
states (stable states lowest and highest in energy, respec¬ 
tively). Further, some of those states can actually be 
unstable, that is no local energy minimum, so that they 
are traversed by coupled jumps. The three simplest de¬ 
fects with the corresponding diffusion mechanisms are 
listed in the following, where for the two charged defects 
obviously also the mirrored case exists: 

— A neutral defect. Its movement is illustrated in 
Fig. a It starts as a nearest-neighbour divacancy, then 
an A atom performs a nearest-neighbour jump into the 
P vacancy, thereby creating a localized A triple defect. 
Successively, the A atom performs a further jump into 
the initial a vacancy, again leading to a divacancy. If 
now a B atom performs the analogous steps, the diva¬ 
cancy has effected a net movement. This is the classical 
divacancy mechanisirPH, but depending on the energet¬ 
ics, one of the triple defect states can be the ground state. 
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FIG. 6. Correlated propagation of doubly charged defect — 
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FIG. 5. Correlated propagation of single-charged defect — 
jump cycle. If the ground states are the single jS vacancy 
states 2 and 8, the illustration depicts the classical six-jump 
cycle linking the ground states 2 and 8, framed by the last 
and first jump of the preceding and subsequent cycles, re¬ 
spectively. If on the other hand the states 1, 3, 7 and 9, each 
with an a vacancy and a P antisite, are the ground states, 
it illustrates the 4-1-2-jump cycle with the introductory two- 
jump defect reorientation, the four-jump cycle, and a further 
subsequent defect reorientation. 

for which case the term triple-defect mechanism has been 
coined^. For implying the opposite case with the divacan¬ 
cies being the ground state, leap-frog mechanism would 
be an appropriate name due to the movement of the two 
constituent vacancies. Independent of these distinctions, 
different jump geometries are possible, which in reality 
will differ in energy. Further, if (say) the A atom jumps 
are much more frequent than B atom jumps because the 
B triple defect is costly in energy, then the defect will 
oscillate rapidly between different orientations of the di¬ 
vacancy state and the A triple-defect state. During this 
time it will shuffle all neighbouring A atoms, before it 
concludes its movement by a B at om jump, thus leading 
to large diffusion asymmetrie d^ l ^^ l 

— A single-charged defect. Its movement by jump cy¬ 
cles is illustrated in Fig. for positive charge, that is, 
starting with a /? vacancy, which in the following will 
be called /? jump cycles. Along a cycle of K nearest- 
neighbour relations, the atoms successively move into the 
vacancy, until K —1 antisite defects are created. The va¬ 
cancy continues on the cycle, so that those anti-structure 
atoms perform a second jump back onto their native 
sublattice. Finally, after 2K — 2 jumps only the initial 
vacancy remains, effectively translated by two nearest- 
neighbour jumps. Also Kj^ atoms and Ar/2 — 1 A 
atoms have been translated by the same amount. For 
the B2 structure, the smallest such cycle correspond to 
K = A, leading to the mechanism being called six-jump 
cycle and the initial estimate that the ratio of the two 
self-diffusion constants is at most 2. In contrast, for in¬ 
stance in the zincblende structure the smallest cycle con¬ 
sists of AT = 6 sites. For the square lattice as illustrated in 
Fig. [^all such smallest cycles are equivalent, whereas for, 
e.g., the B2 structure there exist three geometrically in¬ 
equivalent 4-site cycles. Also here different states can be 
the ground and apex states. Specifically, for significant 
asymmetries the second state with a /3 antisite and an a 
vacancy can be lower in energy than the single /? vacancy. 
In this case the six-jump cycle would be more appropri¬ 


ately called 4-|-2-jump cycle, as the vacancy-antisite pair 
would be able to re-orient freely within its cell, until the 
vacancy eventually performs a cycle during which only 
two additional antisites are created. This had already 
been noted by Athenes and coworker^, but its relevance 
for actual systems as will be discussed in Sect. |IV A] has 
not been recognized yet. 

— A double-charged defect. Its movement for the case 
of an initial a vacancy pair (called a waltzing-step mech¬ 
anism) is illustrated in Fig. It starts with a B atom 
moving into one of the a vacancies. Then an A atom 
performs a double jump via the vacated /3 site into the 
other a vacancy, after which the B atom retakes its ini¬ 
tial site. As a consequence, just as with the percolation 
mechanisms but in contrast to jump cycles and the di¬ 
vacancy mechanism, in this correlated mechanism only 
one constituent experiences a net translation. A fitting 
name for the described correlated movement is waltzing- 
step mechanism, as the B atom stepping temporarily out 
of the way for the A atom to pass perfectly mirrors the 
movement of the female and male partners during one 
bar of the classical ballroom dance. This mechanism has 
never been discussed yet, but below it will be shown that 
under reasonable assumptions it is a priori not less prob¬ 
able than for instance the divacancy mechanism. 


C. Completeness 

Recapitulating the assumptions stated above, viz. a) 
small defect concentrations corresponding to a movement 
of the diffusing complex defect in a perfect lattice and 
b) this movement being composed of nearest-neighbour 
jumps into vacancies, above-discussed three qualitatively 
different mechanisms (divacancy mechanism, jump cy¬ 
cles, and waltzing-step mechanism) constitute an exhaus¬ 
tive list of possible candidates for the diffusion mecha¬ 
nism with lowest activation energy under the additional 
reasonable assumption that c) the introduction of ad¬ 
ditional transient disorder costs additional energy. This 
can be shown as follows: First, all three mechanisms con¬ 
stitute those sequences of jumps where the respective de¬ 
fects or defect complexes can migrate with the smallest 
amount of transient order. This claim can be formally 
verified by surveying the finite number of states that can 
be reached from the respective initial states under the 
constraint to not introduce more antisite disorder than in 
the candidate mechanisms. Further, allowing more than 
two vacancies does not allow the defect complex to mi¬ 
grate with less antisite disorder than both the divacancy 
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and the waltzing-step mechanisms: The movement of a 
defect complex in an otherwise ideal lattice necessitates 
the movement of both constituents, at least a temporary 
stepping aside as in the waltzing-step mechanism. As 
an a site’s direct neighbours are exclusively /3 sites and 
vice versa, any movement of component A by nearest- 
neighbour jumps has to involve an a vacancy and a /? 
antisite, vice versa for B. The divacancy mechanism has 
only a single additional vacancy in the respective corre¬ 
sponding states. Naturally, its apex energy cannot be 
lowered by allowing additional vacancies. If the effective 
defect formation energies should show a large asymme¬ 
try, the waltzing-step mechanism can become active and 
obviate the more expensive transient triple defect, at the 
cost of an additional transient antisite defect. Introduc¬ 
ing a third a vacancy (or even more) to an a waltzing- 
step mechanism does not obviate the transient a anti¬ 
site, while the introduction of a /3 vacancy allows the 
divacancy mechanism to proceed, where the second a 
vacancy can be eliminated. 

Of course, interactions between neighbouring defects 
as well as different energy barriers between the respective 
states can in principle enable mechanisms to be active 
that have not been considered above. Also, mechanisms 
involving interstitial sites (of which direct jumps to va¬ 
cancies farther than on nearest-neighbour sites are the 
most simple representative) cannot be ruled out by such 
considerations. Finally, thermal disorder or significant 
constitutional defect concentrations below the percola¬ 
tion threshold will lead to a blurring of above clear-cut 
distinctions, where diffusion over small scales is promoted 
by analoga of the anti-structure or vacancy bridge mech¬ 
anisms on non-percolating defect clusters, but a jump 
cycle or a cooperation of two vacancies is necessary for 
the diffusing defect to jump to another cluster. 


D. Estimating energetics 

Here the idea of a given diffusion mechanism’s activa¬ 
tion energy being the sum of non-interacting effective de¬ 
fect formation energies that was used qualitatively above 
will be explored quantitatively. This is the simplest con¬ 
ceivable model that takes into account the point defect 
energetics as described by the non-trivial parameters ^ 
and Tj. Such a model has already been used by Kao 
and Chan^^ for determining whether a or (3 jump cy¬ 
cles dominate in a given system. In addition to equating 
the energy difference between apex and ground state with 
the activation energy as discussed in |III A[ neglecting the 
interaction between neighbouring point defects and as¬ 
suming a constant barrier height separating neighbouring 
stable states, it also considers only jumps of single atoms 
into nearest-neighbour vacancies as opposed to coupled 
jumps of several atoms. However, for the present purpose 
of determining the dominant mechanism this coupling 
will affect the energetics of all mechanisms and therefore 
cancel in the first approximation. 


It is expedient to split the diffusion mechanisms’ ac¬ 
tivation energy into the formation energy of the initial 
state and the increase in potential energy during the 
jump sequences. The former contribution follows from 
summing the respective effective defect formation ener¬ 


gies from Eqs. ([^, resulting in 

=-®o(l /2 + ^/ 2 ) ( 6 a) 

el = Eoil/A + ^/A + 2r^/5 + fi) ( 6 b) 

2el = Eoil/2 + ^/2 + 4r^/5 + 2fi) ( 6 c) 

= Eq{1/2 + 77 / 2 ) ( 6 d) 

+ Me^ = Eq{—^ + rj). ( 6 e) 


for the divacancy, a jump cycles, a waltzing-step, (3 anti¬ 
structure bridge, and /3 vacancy bridge mechanisms, re¬ 
spectively. As the percolating mechanisms are to be con¬ 
sidered only in the pertaining off-stoichiometric cases, the 
corresponding /x has already been chosen so that the ef¬ 
fective formation energy of the M necessary percolating 
defects vanishes, with M a large number proportional to 
the system size. 

There are four elemental transitions the nearest- 
neighour jump mechanisms are made up of, viz. the 
transition of an A atom from its native lattice into a 
j3 vacancy or back, and the analogous transitions for 
a B atom. The energy difference for the transition of 
an A atom from the a to the (3 sublattice is given by 
ea-j ./3 = + e;s - = £^o(l /2 - C /2 + ?7), similarly 

+ e. - £« = ^^o(l/2 - C/2 - v)- These en¬ 
ergy differences will be used for the analysis of the dif¬ 
fusion mechanisms’s intermediate states. Specifically, if 
< Q or r] < —1/2-I-C/2 the generation of a /3 antisite 
defect corresponds to an energy gain, analogously for an 
a antisite defect if 77 > 1/2 — C/2. 

For the divacancy mechanism, the energies of the inter¬ 
mediate states relative to the initial divacancy state are 
(0, 0, 0). As a consequence, if—l/2-|-C/2< 

r] < 1/2 — C/2 the divacancy state is the ground state, 
while a triple-defect state is the ground state if rj is out¬ 
side of this interval. The apex state is always one of the 
two triple-defect state, with an energy difference to the 
initial state of i?o(l/2 — C/2 -I- \r]\). 

The jump cycle will be considered in its form with an 
initial a vacancy for cycles with K = 4 as is appropri¬ 
ate, e.g., for the B2 lattice. The sequence of energies are 
( 0 , 2 ef_^„-ke^^^), followed by a sym¬ 

metrical descent. The second intermediate state is higher 
than the initial state by the positive energy needed to 
form an antisite pair, the same applies for the difference 
of the third to the first intermediate state. Therefore, 
only if ry < 1/2 — C/2 the initial state with the single 
vacancy is the ground state and the third intermediate 
state the apex state so that this mechanism is properly 
called six-jump cycle. In the converse case, the first in¬ 
termediate state is the ground state and the second inter¬ 
mediate state is the apex state, so that diffusion by this 
mechanisms corresponds to sequences of two-jump reori- 




10 



77 <- 1/2 +C /2 

-1/2 +C/2 < 77 < 1/2-C/2 

1/2 - C/2 < 77 

divacancy 
a jump cycles 
a waltzing-step 
/3 anti-structure bridge 
P vacancy bridge 

I + I 77 I 

7/4 - 5C/4 - 377/5 + /i 

1 - 77/5 + 2p 

C/2 - 77/2 

I + I 77 I 

7/4 - 5C/4 - 377/5 + + 

3/2 - C/2 + 477/5 + 2p 

1/2 + 77/2 
1/2-3C/2 + 277 

1 + I 77 I 

5/4 - 3C/4 + 277/5 + + 
3/2 - C/2 + 477/5 + 2p 
1/2 + 77/2 
1/2-3C/2 + 277 


TABLE IV. Activation energies (effective formation energy of initial state plus increase towards apex state) of the considered 
diffusion mechanisms in units of Eq. The expressions for the mirror cases {a o /3) correspond to swapping the signs of 77 and 


entations of the vacancy-antisite pair alternating with a 
four-jump cycle, which is more fittingly called 4-|-2-jump 
cycle, as discussed already above. 

In the case of the waltzing-step mechanism with two 
initial a vacancies, the sequence of energies up to the 
state with mirror symmetry is 

Either the first or the second intermediate state are the 
apex state depending on whether rj < —1/2 -1-^/2 or not. 

The anti-structure bridge mechanism under A-excess 
alternatingly passes through states with energy 0 and 
relative to the initial state with an a vacancy, 
and either of those can be ground and apex states. For 
the vacancy bridge mechanism the energy difference of 
the intermediate state is which is positive in the 

region of constitutional vacancies. 

The apex state energies of the respective mechanisms 
and therefore their activation energies are tabulated in 
Tab. |IV[ The resulting dominant mechanisms defined by 
the smallest activation energies are illustrated in Fig. 

As in the case of point defects, also here all of the mech¬ 
anisms considered above are active for some (C, ??), and 
furthermore for each mechanism there are representatives 
among the considered systems where it is expected to be 
active. The critical points are (3/7, ±1/7) and (3/5,0) 
for stoichiometry, (0,0), (1/2,0) and (5/7,1/7) for B 
diffusion under A-excess, (1/2,0) and (1/3,1/3) for A 
diffusion, and (—1,—1/3), (—1/3,0) and (0,1/2) for the 
vacancy bridge mechanism. Verifying the presented as¬ 
signments of mechanisms to regions in (^, ? 7 )-space is con¬ 
ceptually simple but tedious and will not be performed 
here. 

In the stoichiometric case, jump cycles and waltzing- 
step mechanisms share the minimal activation energy 
over a large region of the phase space, for instance (3 
jump cycles with the a waltzing-step mechanism in the 
A triple-defect region according to Fig. with 77 < 0. 
This is because the apex states of the respective mech¬ 
anisms agree up to a single point defect, which is an a 
vacancy in the case of the a waltzing-step mechanism 
and a /3 antisite for f3 jump cycles, having equal effective 
formation energies in the case of A triple defects. The 
waltzing-step mechanism, which leads only to diffusion 
of one component, is never exclusively active, so that 
diffusion of both constituents will be coupled. A point 
worth noting is that the major part of this region fulfills 
\'q\ > 1/2 — ^/2, so that the corresponding jump cycles 


correspond to 4±2-jumps, with the ground state made up 
of a q; vacancy and an /? antisite. On the other hand, the 
divacancy mechanism is active only for \r]\ < 1/2 — ^/2, 
where it will take the form of the vacancy-leapfrog mech¬ 
anism. Its dominance in the region of negative ^ and 
small absolute values of rj can be easily understood, as 
small ^ favour vacancies compared to the other mecha¬ 
nisms, which rely more on antisites, while 77 shifts barrier 
height between the two triple-defect states, so that the 
smallest apex energy results from 77 = 0. 

In the off-stoichiometric case, the degeneracies between 
jump cycles and the waltzing-step mechanisms are lifted. 
Specihcally for A-excess, for 77 < 0 the /3 jump cycles 
dominate, while for 77 > 0 the /3 waltzing-step mechanism 
dominates, causing diffusion of only the B minority com¬ 
ponent. This leads to the result that, with the exception 
of FeCo, all systems considered here show coupled dif¬ 
fusion for both constituents in the transition metal-rich 
case, while in the converse case the minority component 
(i.e., the transition metal) should diffuse with a smaller 
activation energy and therefore faster. 

It can be seen that for percolating constitutional de¬ 
fects, the anti-structure bridge mechanism is uncondi¬ 
tionally responsible for majority diffusion, while the two 
antisite defects in the apex state of the vacancy bridge 
mechanism limit its effectiveness to a strip along the main 
diagonal. Still, according to the present model for all 
considered systems except for In-rich Pdin percolating 
mechanisms will dominate majority diffusion once they 
are made possible by sufficient deviations from stoichiom¬ 
etry. 

There exists a region in the phase space where a perco¬ 
lation mechanism dominates majority diffusion, while the 
waltzing-step mechanism is responsible for minority dif¬ 
fusion. In this case, the information plotted in Fig. [^can¬ 
not decide which component will have a smaller activa¬ 
tion energy of diffusion. Pertinent calculations show that 
for the case of constitutional antisites, the anti-structure 
bridge mechanism corresponds always to the lower ac¬ 
tivation energy, while the vacancy bridge mechanism is 
lower in energy than the waltzing-step mechanism if and 
only if 77 < 1/4 ± ^/4. 

While above paragraphs discuss which component dif¬ 
fuses faster for given system and composition, the com¬ 
plementary question as to how the activation energy 
varies with composition is answered by Fig. It shows 
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FIG. 7. Dominant diffusion mechanisms for non-interacting defect energetics. The divacancy mechanism is orange, jump cycles 
yellow (a) and green (/?), waltzing-step mechanisms pink (a) and blue (/?). Striped areas denote degeneracies. The dashed 
lines correspond to \ri\ = 1/2 — ^/2, which decides on the apex states as discussed in the text. The left-top panel treats 
the stoichiometric case, the other three panels A-rich off-stoichiometry, namely B diffusion (right top), A diffusion without 
percolation of constitutional defects (left bottom), and A diffusion including percolating mechanisms (right bottom), where 
dark gray corresponds to the anti-structure bridge mechanism and light gray to the vacancy bridge mechanism. 


that for major parts of the phase space the activation 
energy of diffusion for a given constituent is smaller if it 
is the minority species. Naturally, in systems where the 
divacancy mechanism is dominant, there is no effect of 
composition on the activation energy, as a divacancy is a 
neutral defect and its concentration is therefore indepen¬ 


dent of /r. Only in a small region of the phase space a 
component’s activation energy of diffusion drops with an 
increase in its concentration. However, once the specific 
percolation threshold is reached, the concerned compo¬ 
nent in most cases can diffuse with the smallest activation 
energy. The additional critical points are (1/7, 3/7) and 
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FIG. 8. Behaviour of activation energy of A diffusion as func¬ 
tion of composition. Yellow denotes smaller activation energy 
for A-excess than for A-deficiency, blue the converse case, and 
orange equal activation energies for both cases. Active per¬ 
colation mechanisms would enlarge the yellow region towards 
the dashed line. 


(-1/13,5/13). 

The proposed model of estimating activation energies 
for the different mechanisms is the simplest conceivable 
model that covers all mechanisms composed of atomic 
jumps into nearest-neighbour vacant sites. In reality, al¬ 
ready in the small-temperature limit effects not treated 
here will complicate the picture. These include inter¬ 
actions between the point defects introduced during the 
various steps of the mechanisms as well as jump barri¬ 
ers separating the stable states that depend on the envi¬ 
ronment. As a consequence, for instance the three geo¬ 
metrically inequivalent six-jump cycles will have different 
saddle point energies (see the discussion in Sect. IVA|. 
Further, in order to avoid especially costly intermediate 
states such as an aluminium vacancy in the case B2 tran¬ 
sition metal-aluminides, the system can perform coupled 
jumps of several atoms. However, this would affect all 
mechanisms, so that for the present purpose of compar¬ 
ing the activation energies of different mechanisms it can 
be expected to cancel to a large amount. Finally, the 
interplay of different mechanisms due to equilibrium dis¬ 
order at high temperature can in practice only be treated 
by Monte Carlo simulations (e.g. Ref. III). In contrast to 
such simulations where a detailed understanding of the 
observed effects can be hard or even impossible to ob¬ 
tain, the reasoning presented here is rather transparent. 
In spite of its shortcomings as discussed above, the qual¬ 
itative conclusions can be expected to be sound. Specifi¬ 
cally, it motivates the existence of the jump mechanisms 


discussed and introduced in Sect. niBi 


E. Pair interactions 

In the model investigated in Sects. [ir^and lllT^ the 
energies of specific states depend only on the occupations 
of the distinct sublattices. In an interaction picture, it is 
the limit of a mean-field model for small defect concen¬ 
trations and therefore corresponds to an infinite interac¬ 
tion range, distinguishing only between sublattices. In 
previous investigations (e.g., Ref. [8| often the opposite 
extreme of nearest-neighbour pair interactions has been 
considered, where the energies follow from summing the 
occupation-dependent contributions exY of (undirected) 
nearest-neighbour pairs and the distinction of sublattices 
arises only due to spontaneous symmetry breaking of the 
thermodynamic ground state. Equating the formation 
energies of isolated point defects leads to the relations 

CAB = Eoi-1 - i)/2Z (7a) 

eAA = Ao(-e + 77)/^ (7b) 

eBB = Ao(-C - rj)lZ, (7c) 

where Z is the coordination number. This implies that 
a pair-interaction description as conventionally used in 
Monte Carlo simulations has a priori the same generality 
as the model used here. However, allowing only a single 
vacancy in the simulation as is often don^ corresponds to 
neglecting one degree of freedom and to potentially sup¬ 
pressing diffusion mechanism relying on the correlated 
movement of several vacancies, i.e., the divacancy and 
waltzing-step mechanisms. 

Different from the model of non-interacting defect en¬ 
ergetics, nearest-neighbour pair energetics lead to effec¬ 
tive interactions between neighbouring point defects 


v„-v^ 

: 2eAB=Eo{-l-0/Z 

(8a) 

Bq-A^ 

■ 2eAB — gaa — sbb = Eo{—l + i)/Z 

(8b) 


'■ Gab ~ gaa = Ao(—I + C ~ 2vj)l2Z 

(8c) 


'■ CAB — cbb = Ao(—1 + C + 2??)/2Z, 

(8d) 

which are 

typically of attractive nature. These 

expres- 


sions specify how the energies of the diffusion mecha¬ 
nisms’ intermediate states have to be modified for per¬ 
forming an analysis as done in Sect. |HID[ but here for 
the assumption of nearest-neighbour pair energetics. 

The resulting phase diagrams are given in Fig. As 
expected, the situation is more complicated than for non¬ 
interacting defect energetics, which is mostly a conse¬ 
quence of the cross-over from one apex state to the next 
not coinciding for the respective mechanisms due to the 
different impact of point defect interactions. The illus¬ 
trated topology is representative for Z > 3, i.e., for all 
relevant cases. Qualitatively, the introduction of local 
interactions (corresponding to Z < oo) leads to a lift¬ 
ing of the degeneracies and the appearance of additional 
phase regions in the outer regions of the phase space. In 
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FIG. 9. Dominant diffusion mechanisms for pair-interaction energetics with Z — 8 a,t stoichiometry (left) and A-excess (right). 
Colour code is the same as in Fig. Dashed lines delineate regions neglecting the waltzing-step mechanism, i.e., indicate the 
mechanism responsible for diffusion of the slower component. 


the inner regions, the phase boundaries do not change 
much, however, so that the non-interacting classification 
(Fig. and the pair-interaction classification (Fig. 
agree on all considered B2 systems with the exception 
of A1 diffusion in Al-rich PdAl and CoAl, which should 
proceed by the divacancy mechanism according to the 
non-interacting model, and by the six-jump cycle start¬ 
ing from the transition metal vacancy according to the 
pair-interaction model. Note that with local interactions 
neighbouring point defects can affect the energetics even 
if they do not participate in the jump, so that the perco¬ 
lating mechanisms cannot be assigned a definite activa¬ 
tion energy. 

The fact that the dominant mechanisms as predicted 
by the opposite extremes of nearest-neighbour and infi¬ 
nite interaction range energetics agree so well gives sup¬ 
port to the models. It shows that interactions between 
neighbouring point defect affect the energies of the inter¬ 
mediate states approximately equally for all mechanisms, 
so that these modifications cancel to a large extent. 


IV. APPLICATION TO SELECTED SYSTEMS 

Here the predictions on the active diffusion mecha¬ 
nisms according to above models will be assessed by com¬ 
paring to calculations and experimental diffusion data on 
NiAl and CoGa, two prototypical and well-investigated 
B2 systems. 


A. NiAl 

Both the non-interacting (Fig. and the pair- 
interaction model (Fig. 1^ predict A1 4-|-2-jump cycles to 
be the dominant mechanism of diffusion in all cases with 
the exception of Ni-deficiency, where the Ni six-jump cy¬ 
cles have comparable energy and the double Ni vacancy 
waltzing-step mechanism should allow faster Al-diffusion. 
Once the respective thresholds are reached, Ni-diffusion 
should proceed by the anti-structure bridge mechanism 
under large Ni-excess and Al-diffusion by the vacancy 
bridge mechanism under large Ni-deficiency. 

For NiAl there have been a number of calculations of 
diffusion energetics. For consi stency, here the results re¬ 
ported by Xu and Van der VerP^^EHl obtained by nudged 
elastic band density-functional theory calculations will be 
reported and used for deriving estimates for mechanisms 
that have not been considered yet. A consistent feature 
is the appearance of coupled jumps involving A1 atoms. 
An A1 vacancy is energetically expensive, therefore it is 
favourable that the vacated A1 site is immediately reoc¬ 
cupied by another atom. If the effective translation of the 
vacancy is along (100), the angle between the respective 
atomic translation vectors is acute in the idealized pic¬ 
ture. Due to short-range repulsion between the jumping 
atoms, the A1 site can therefore not be covered as effec¬ 
tively as for (110) jumps with an obtuse angle or even 
collinear (111) jumps. This is likely the reason for the 
observed increase of barrier energies for equivalent cou¬ 
pled jumps from (111) over (110) to (100) as observed 
already in Ref. 01] and can be assumed to be of general 

















14 


validity. 

For the Ni vacancy six-jump cycle Xu and Van der 
Ven give an apex energy of 2.59 eV above the single va¬ 
cancy state for the (110) configuration, and higher values 
for the (100) configurations or the equivalent mechanism 
with a second-nearest neighbour A1 jump through the 
apex. Along with all other such calculationPiH^ the A1 
vacancy six-jump cycle (or rather 4-|-2-jump cycle as will 
be seen) is not considered, which is explained either by 
its claimed notably higher activation energie^^ or not 
justified at all. However, an estimate of its apex en¬ 
ergy can be put together from those of its components: 
The nearest-neighbour pair of a Ni vacancy and a Ni 
anti-structure atom is 0.6 eV lower in energy than an iso¬ 
lated A1 vacancy, and the apex energy of this transition 
is 0.35 eV with respect to the initial vacancy. There¬ 
fore the vacancy-antisite pair is the ground state, and 
it can perform its two-jump reorientations with 0.95 eV. 
The four-jump part of the mechanism then proceeds by 
a coupled (110) Al-Ni-jump into the Ni vacancy, which 
is just the first double jump of the (110) Ni vacancy six- 
jump cycle with 2.36 eV. However, the presence of the 
Ni anti-structure atom can be expected to decrease this 
energy further, just as it decreases the barrier energy of 
the coupled Al-Al-jump from 1.01 eV in the isolated case 
to 0.53 eV in the Ni six-jump cycle. Therefore, the acti¬ 
vation energy of the A1 4-|-2-jump cycle can be expected 
to be 1.76 eV above the isolated A1 vacancy, probably 
decreased by a few tenths of eV. Note that while the 
four-jump part as described here leads to a (110) dis¬ 
placement of the vacancy, in the six-jump description it 
would correspond to a (100) cycle. 

The divacancy mechanism should take the triple-defect 
modification with the Ni triple defect as ground state. 
For the (100) jump its apex energy during the A1 jump is 
0.75 eV above the nearest-neighbour divacancy, the en¬ 
ergy of which should not deviate by more than 0.1 eV 
from isolated vacancieJi^, while the energy is much higher 
for other jump geometries. For the double Ni vacancy 
waltzing-step mechanism the energetics have not been 
computed yet, too. Again the coupled (110) Al-Ni-jump 
with 2.36 eV can be considered as an estimate for its en¬ 
ergy with respect to the double Ni vacancy, with perhaps 
slightly smaller energy for (111) jumps but 2.88eV for 
(100) jumps. The presence of the second Ni vacancy will 
probably affect these barrier energies to a small degree. 
The converse case of an A1 waltzing step corresponds to 
the 0.75 eV of an A1 jump into a divacancy relative to 
the effective formation energy of the complex of a di¬ 
vacancy and a Ni anti-structure atom. To compare, the 
direct second-nearest neighbour jumps within one sublat¬ 
tice cost 2.76 eV for Ni and 1.49 eV for Al, which shows 
that the availability of vacancies with low or vanishing 
effective formation energies decreases the energy for di¬ 
rect jumps within a sublattice by a non-jumping atom 
withdrawing into the vacancy, which is the essence of the 
waltzing-step mechanism. 

The barrier energy for the Ni anti-structure bridge 


mechanism is around 1 eV relative to the Ni vacancy, 
depending on the surrounding configuration of the con¬ 
stitutional antisites. The Al vacancy bridge mechanism 
proceeds via coupled jumps with 1.65 eV for (100) an¬ 
tisite translations and 1.01 eV for (110), again probably 
even less for (111). Constitutional defects show repul¬ 
sions of about 0.1 eV over (100), while interactions over 
farther distances are negligible^S. With the resulting sup¬ 
pression of (100) coordinations the probabilities for (110) 
and (111) pairs will rise, so that the respective percola¬ 
tion thresholds of cni = 0.55 and 0.475 will probably be 
even smaller. 


As Xu and Van der Ven do not report their defect 
formation energies, an average of the values in Tab. |HI| 
with Eq = 4.14 eV, ^ = 0.228 and rj = —0.460 will be 
assumed in the following. Computing the activation en¬ 
ergies of the various mechanisms by adding the corre¬ 
sponding effective formation energies for the necessary 
vacancies to the barrier energies gives the following hi¬ 
erarchy: For Ni-deficiency, the Al vacancy bridge mech¬ 
anism has an activation energy of 1.97 eV, followed by 
the Ni waltzing-step mechanism with about 2.36 eV, the 
Ni six-jump cycle with 2.59 eV and direct Ni jumps to 
second-nearest neighbours with 2.76 eV. Further options 
like the triple-defect mechanism have 3.29 eV or higher 
and can therefore safely be neglected. For Ni-excess, the 
Ni anti-structure bridge mechanism is lowest in energy 
with 2.11 eV, followed by direct Al second-nearest neigh¬ 
bour jumps with 2.91 eV and the Al 4-|-2-jump cycles 
with a few tenths of eV below 3.18 eV. The triple-defect 
and Al waltzing-step mechanisms with 3.29 eV are al¬ 
ready clearly higher. For stoichiometry, the energies of 
Ni six-jump cycles, Al 4-|-2-jump cycles, the triple-defect 
mechanism, and direct Al jumps all cluster closely around 
3.3 eV. These values reproduce the rankings of the sim¬ 
ple models of Sects. HID and HI E[ only the dominance 
of the Ni six-jump cycle for Al diffusion at Ni-deficiency 
due to the inexpensive coupled Al-Al-jumps is missed, 
and the second-nearest neighbour jumps are obviously 
not included. 


The experimental evidence in this system convincingly 
demonstrates constant Ni tracer diffusivities as a func¬ 
tion of composition for Ni-deficiency and smoothly rising 
diffusivities with Ni-excess. This is mirrored in the ac¬ 
tivation energies of Ni diffusion with a constant value of 
3.0 eV on the Ni-deficient side and sharply dropping over 
about 53 % Ni content!^. Assuming In diffusion to mir¬ 
ror Al tracer diffusioiP^, it shows a sharp minimum in 
diffusivity and a maximum in activation energy around 
stoichiometry. In contrast, the simulated diffusivities for 
both component^ show a sharp minimum at stoichiom¬ 
etry, a steep increase towards Ni-excess, and a smaller 
increase with equal diffusivities of Ni and Al towards Ni- 
deficiency. 

A number of open questions remain for diffusivity in 
this system: First, the concentration of the bound di¬ 
vacancy, being a composition-conserving excitation, is 
nearly independent of compositiorP^. This will also hold 
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for the migration energy of the A1 jump into the diva¬ 
cancy, therefore the contribution of the triple-defect to 
diffusion is independent of composition. The decom¬ 
position of diffusivities in Ref. 0 shows that when go¬ 
ing from stoichiometry to Ni-deficiency, as expected, the 
A1 vacancy bridge jumps followed by Ni six-jump cycles 
and the Ni waltzing-step (which remained unrecognized 
and probably classihed as coupled Ni-Al back and forth 
jumps) overtake the triple-defect mechanism and there¬ 
fore yield a higher contribution to diffusion in absolute 
values. As the composition is in the region where the va¬ 
cancy bridge mechanism should be able to contribute to 
long-range diffusion, it is surprising that specihcally the 
simulated A1 diffusivity does not grow much compared 
to the stoichiometric case, and also the Ni diffusivity is 
comparable to the values at small Ni-excess. For larger 
Ni-excess the A1 4-|-2-jump cycles should become active, 
with its two-jump reorientation classified as Ni back and 
forth hops and the four-jump cycle again as coupled Ni- 
Al-hops. However, the expected contribution from direct 
A1 second-nearest neighbour jumps remain unaccounted 
for. 

On the other hand, the experimental finding of 
smoothly rising Ni and A1 diffusivities with Ni-excess 
can be explained by the Ni anti-structure bridge jumps, 
direct A1 jumps, and A1 4-|-2-jump cycles, all of which 
will rise in frequency with Ni-excess due to entropic rea¬ 
sons. Spe cifical ly the latter mechanism has hitherto been 
overlookecP^I^. While in the Ni-deficient case the rising 
A1 diffusivity will probably be due to the A1 vacancy 
bridge mechanism, either solely responsible or enhancing 
the A1 diffusivity in Ni six-jump cycles, the apparently 
constant Ni diffusivity and activation energy is perhaps 
only an accidental effect and can arise from a gradual 
shifting from an anti-structure bridge assisted A1 4-1-2- 
jump cycle to a combination of Ni waltzing-step mech¬ 
anism, Ni six-jump cycles and direct Ni second-nearest 
neighbour jumps, as due to high-temperature disorder 
the chemical potential /r and therefore the point defect 
concentrations show a smooth behaviouJ^. This seems 
plausible as all of these three mechanisms are expected 
to have a lower activation energy than the triple- defect 
mechanism assumed to be responsible previouslj^^sEll^ 
which will more than balance their entropic disadvan¬ 
tage due to the reliance on constitutional vacancies. For 
the relation between the high-temperature experimental 
activation energies and zero-temperature calculations see 
also the following section. 


B. CoGa 

This compound is the prototypical example for a sys¬ 
tem where diffusion is assumed to proceed via the triple¬ 
defect mechanism. In tracer diffusion measurements over 
a large temperature range, deviations from linearity in 
the Arrhenius plots of both constituents’ diffusivities 
were observed^. It turned out that the experimental data 


could be well fitted assuming two or more activation ener¬ 
gies per component and composition, ranging from 2 eV 
to 5 eV, with the additional restriction that the larger Co 
activation energy be equal to the smaller Ga activation 
energy. This was taken as evidence for the simultaneous 
operation of several distinct diffusion mechanisms, where 
the one corresponding to the coupled activation ener¬ 
gies was postulated to be the triple-defect mechanisnP 
in analogy to the case of pure elements, where a contri¬ 
bution from divacancies was assumed to be responsible 
for deviations from linear Arrhenius plots. 

However, also the contrasting opinion that anhar- 
monicity of a single kind of defect can be responsible 
for the observed curvatures has long been helcP^. Only 
recently the analogous problem of curvatures in the Ar¬ 
rhenius plots of vacancy concentrations has been settled 
by ab-initio calculations taking anharmonicity fully into 
accounlP^, where it has been shown that the formation 
entropy in cases as simple as single A1 and Cu vacancies 
depends signihcantly on temperature, being able to fully 
explain the experimentally observed curvatures. 

In the case of an ordered compound, a temperature- 
dependent order parameter can be a further source of 
curvature.^ In the light of this, the initial justification 
for the triple-defect mechanism to be dominant at specihc 
temperature ranges in CoGa seems doubtable, given also 
the fact that the respective diffusivity prefactors I?o of 
the multi-exponential fits span ten orders of magnitude. 
Indeed, fitting the diffusion dataP by a single mechanism 
that has an entropy varying linearly with temperature 
leads to values of per degree of freedom that are prac¬ 
tically equal to those of the multi-exponential fits. For 
the Co-rich and stoichiometric compositions, both con¬ 
stituents’ diffusion energies range from 2.0 eV to 2.7 eV, 
while the fitted values become physically implausible at 
Al-rich compositions. This shows that at least for stoi¬ 
chiometric and Co-rich compositions the observed curva¬ 
ture cannot be taken as indication of a crossover between 
two diffusion mechanisms. 

In default of point defect energy calculations it seems 
safe to assume that the situation here corresponds most 
closely to NiAl, also displaying constitutional transition- 
metal vacancies albeit with a smaller Eg leading to higher 
disorder at comparable temperature^^ and a wide region 
of B2 phase stability. Furthermore, also the experimental 
behaviour of the diffusivities conforms, with Ga showing 
faster diffusion on either side of stoichiometry, while the 
Co diffusivity seems unaffected by Ga-excess. The dif¬ 
fusion acceleration due to deviations from stoichiometry 
is smaller than in the case of NiAl, probably because of 
the higher degree of disorder. As in NiAl, contrary to 
expectations there is no discernible setting-in of the anti¬ 
structure bridge mechanism at some Co-rich percolation 
threshold, rather this mechanisms seems only to play a 
secondary role in enhancing the Co diffusivity effected 
by some other mechanism, active on both Co and Ga, 
with progressive Co-excess. The reason for this proba¬ 
bly again lies in the comparatively low degree of order, 
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which agrees with Co diffusion as measured by quasi¬ 
elastic neutron scattering being well described by uncor¬ 
related jumps between both sublattice^^. On the other 
hand, the fact that Ga diffuses decidedly faster than Co 
in Co 45 . 2 Ga 54.8 at low temperatures is an even clearer in¬ 
dication for the vacancy bridge mechanism than in NiAl, 
which is probably the most relevant feature displayed by 
the available diffusion data in this compound. 

V. CONCLUSIONS 

To conclude, in this article a complete characterization 
of the defect thermodynamics in binary systems with two 
chemically inequivalent sites has been given. The various 
possible combinations of constitutional defects and ther¬ 
mal excitations have been enumerated, and specifically 
the relation of constitutional defects and thermal triple¬ 
defect excitations has been clarified. Published compu¬ 
tations of point defect energetics have been classified ac¬ 
cordingly, and it has been demonstrated that all the the¬ 
oretical possibilities have representatives among actual 
binary compounds. 

Further, a list of jump mechanisms composed of 
nearest-neighbour translations of atoms into vacancies 
has been given, and its exhaustiveness for small defect 
concentrations has been motivated. These mechanisms 
are the well-known six-jump cycles and divacancy mech¬ 
anism, while in addition the waltzing-step mechanism 
has been proposed here for the first time. For percolat¬ 
ing clusters of antisites and vacancies, the anti-structure 
bridge and the as yet little known vacancy bridge mecha¬ 
nisms, respectively, have been quoted, along with precise 
values for the appropriate percolation thresholds. 

Two models for computing diffusion energetics directly 


from point defect energies have been considered, one in 
terms of non-interacting defects and the other for nearest- 
neighbour pair interactions. As for the point defect ther¬ 
modynamics, a classification according to dominating dif¬ 
fusion mechanisms at low temperatures has been given 
for both models. The close agreement between their pre¬ 
dictions is judged as a strong argument for their quali¬ 
tative correctness. Specific findings are that the hitherto 
not considered waltzing-step mechanism will dominate 
minority diffusion over a large range of systems, and that 
the six-jump cycles will in most cases actually proceed 
as 4-|-2-jump cycles. This has the consequence that for 
transition metal-group III compounds such as NiAl the 
six-jump cycles starting from the main group vacancies 
are expected to dominate, contrary to prevailing opinion. 
Also, the divacancy mechanism is expected to be less rel¬ 
evant for strongly asymmetric systems, where it would 
take the triple-defect form, than for symmetric systems. 

Finally, the model predictions have been compared 
to computed diffusion energetics in NiAl, and satisfac¬ 
tory agreement has been demonstrated. These findings 
lead to new interpretations of experimentally determined 
composition-dependent diffusivities in NiAl and CoGa. 
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